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This theoretical paper investigates the formation of ground state molecules from ultracold cesium 
atoms in a two-color scheme. Following previous work on photoassociation with chirped picosecond 
pulses [Luc-Koenig et al., Phys. Rev. A 70, 033414 (2004)], we investigate stabilization by a second 
(dump) pulse. By appropriately choosing the dump pulse parameters and time delay with respect to 
the photoassociation pulse, we show that a large number of deeply bound molecules are created in the 
ground triplet state. We discuss (i) broad-bandwidth dump pulses which maximize the probability 
to form molecules while creating a broad vibrational distribution as well as (ii) narrow-bandwidth 
pulses populating a single vibrational ground state level, bound by 113 cm~^. The use of chirped 
pulses makes the two-color scheme robust, simple and efficient. 



Oh: 



O 
Oh' 



> 

o 

o^ ; 
o . 
oo 

o 

o : 
o 

>^ 

^ ; 

Oh- 



X 



PACS numbers: 33.80.Ps,33.80.-b,32.80.Qk,33.90.+h 



I. INTRODUCTION 

The field of cold molecules is rapidly developing, since 
many research groups are working at the production of 
dense samples of cold (T <1 K) or ultracold (T < 100 
fiK) molecules ,lj. Moreover, the field of molecular con- 
densates has opened recently, with experimental evidence 
for molecules formed in a degenerate gas, starting either 
from an atomic Bose-Einstein condensa te B , B , B, [3, B , 
or from an atomic Fermi gas 7,, ^, ^ ilOl IllL Il2 . uM- 
This intense activity is motivated by possible applica- 
tions to very precise measurements (lifetimes, checks on 
parity violation, determination of a possible dipole mo- 
ment of the electron as a check for the standard model 
in elementary particles), Bose-Einstein condensation of 
complex systems (molecule laser), or the emergence of a 
new ultracold chemistry. 

Among the various routes to produce ultracold 
molecules, photoassociation has emerged as a purely op- 
tical technique starting from an assembly of laser-cooled 
atoms in a trap. The reaction takes place when a pair of 
ground state atoms interacting via the potential V^iound 
absorbs a photon red-detuned compared to the atomic 
line. At resonance, one vibrational level in an excited 
electronic potential Kxc of the molecule is populated. 
The reaction can be interpreted as a vertical transition at 
large distances, forming a long range molecule. We shall 
consider the example of two cold cesium atoms colliding 
in the lower triplet a^E^(6s -I- 6s) potential and forming 
a molecule in a vibrational level v' of the 0~(6s -I- 6^3/2) 
potential. The first experimental observation of ultracold 
Cs2 molecules has made use of this reaction [l4] . 

In order to make molecules in the ground or lower 
triplet state, the photoassociation step must be followed 
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by a stabilization step. The photoassociated molecules 
are short-lived and decay via spontaneous emission, most 
often through a vertical transition at large interatomic 
distances R giving back a pair of atoms. However, spe- 
cific mechanisms which favor radiative decay at shorter 
distances exist and allow for populating bound vibra- 
tional levels in the potential Vground(^) (or possibly, for 
instance in case of heteronuclear dimers, to another po- 
tential correlated to the same asymptote 15]). In the 
special case of the double well 0j{6s + 6^3/2) excited po- 
tential of Cs2, the smooth barrier between the two wells 
acts as a "speed bump" and slows down the vibrational 
motion in the external well such that vertical transitions 
can take place at intermediate distances i?int- Stable 
molecules formed by spontaneous emission through this 
mechanism have been detected in a MOT ■ Other effi- 
cient stabilization mechanisms such as resonant coupling 
111 111 El exist. 

In a condensate, in order to retain coherence, sponta- 
neous emission should be avoided and a stabilization step 
involving stimulated emission should be implemented 
within a two-color experiment. However, with continu- 
ous wave (cw) lasers such experiments are hardly feasible 
since the scheme is fully reversible: Population is pumped 
back from bound levels of the ground state to the excited 
state, and then to continuum levels of the ground state, 
i.e. the molecule is dissociated. The use of laser pulses 
becomes thus unavoidable. 

In previous work 0, I20I l2ll ]. the possibility to in- 
crease the efficiency of the photoassociation reaction by 
use of chirped laser pulses was considered. The frequency 
of such pulses has a linear time-dependence. The pa- 
rameters of the pulse were optimized to achieve a total 
transfer of population under adiabatic following condi- 
tions m 111 III. It was shown HO, 113 that a significant 
improvement of the photoassociation rate compared to 
cw excitation can be obtained, since during the pulse the 
resonance condition is fulfilled for several vibrational lev- 
els v' . However, the quantity which must be improved is 
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the formation rate of ground state molecules rather than 
the photoassociation rate. One advantage of excitation 
with chirped pulses is that the chirp parameter can be 
chosen to modify the shape of the vibrational wave packet 
formed in the potential Vcxc- This may improve the ef- 
ficiency of the stabilization step, and thus the formation 
rate of ground state molecules. 

The aim of the present work is therefore to suggest a 
two-color experiment using short laser pulses and to in- 
vestigate how efficiently a second pulse can transfer the 
population to the bound vibrational levels of the ground 
state. Starting from the vibrational wave packet formed 
in the excited state by photoassociation, we introduce a 
second pulse which is delayed in time from the first pulse. 
The time delay is chosen such that the wave packet has 
reached the vicinity of the inner turning point, where a 
vertical radiative transition to the ground state is ex- 
pected to be efficient. The parameters of the second 
pulse, as well as the time delay between the two pulses are 
varied to improve the efficiency of the transfer, with two 
different goals: to either maximize the total population 
transfer, or to selectively populate one vibrational level. 
The paper is organized as follows: In Sec. ^l^e recall 
the main results of previous work on photoassociation 
with chirped laser pulses. The two-color scheme is intro- 
duced in Sec. IIIII and the best choice for the dump pulse 
discussed by considering the Franck-Condon overlap be- 
tween the vibrational wave packet in the excited state 
and wave functions of the bound levels in the ground 
state. The results for maximum population transfer are 
presented in Sec. II VI while Sec. ^reports on the selective 
population of a single level. Finally Sec. I VII concludes. 



II. CHIRPED PULSE PHOTOASSOCIATION IN 
CS2 AND CREATION OF A FOCUSSED 
VIBRATIONAL WAVE PACKET 

We briefly summarize the main results of Refs. [20II2H 
for the photoassociation of cold cesium atoms by chirped 
laser pulses. The main difference to cw excitation lies 
in two facts: Several levels can be populated in the 
photoassociation step, and the vibrational wave packet 
can be shaped in order to allow for spatial "focussing", 
as described below. We start from a transform-limited 
pulse with Gaussian envelope, hereafter referred to as 
pump pulse, with maximum at t = tp, and with full 
width at half maximum of the temporal intensity profile 
(FWHMI) tl — 15 ps. Chirping the pulse stretches the 
FWHMI to rc= 34. 8 ps and decreases the maximum in- 
tensity II ~> II-^ {II =120 kW cm~^), conserving the 
energy carried by the field. Moreover, 98% of this energy 
is carried during the time window [tp — Tc,tp + tc]. We 
consider the reaction 

2Gsi6s, F^4:) + h[LJL + Xp{t-tp)]^ (1) 
Cs2(0;(6s 2^1/2 + 6p ^P3/2;v', J' = 0)) 



where the carrier frequency is chosen red-detuned 
relative to the D2 atomic line, at resonance with the 
v'=v']^—98 level in the external well of T^xc which is bound 
by 2.65 cm^^. xp is the linear chirp rate in time domain, 
chosen negative and equal to -0.025 cm~^ ps~^. The pa- 
rameters of the pump pulse, labelled V- are recalled in 
Table m We will also consider other possibilities for the 
pump pulse, with positive or zero chirp, or with a smaller 
detuning: such pulses are respectively labelled V+, Vo 
and V]^^. The latter pulse has a central frequency reso- 
nant with the excited level w'=122. The parameters for 
V}^^ have been chosen to ensure the focussing condition 
without populating continuum levels in the photoasso- 
ciation step. Note that in order to compare to earlier 
papers, we have also considered another the temporal 
width, 

tS = V2tc, 
fp = rc/TL = T<?/r^, (2) 

where is now defined as the full width at half 
maximum of the electric field profile (FWHME) of the 
transform-limited pulse, while in Refs. Ij^, 0| tl was 
defined with respect to intensity (FWHMI). fp is the 
stretching factor due to chirping. During the time win- 
dow, the central frequency is swept from ujl — X^c to 
w^-l-XTc such that the resonance condition is satisfied for 
15 vibrational levels around t;^=98 in the energy range 
2^ I X I Tc = 2 X 0.87 cm~^. Depending upon the sign 
of the chirp, the frequency is increasing {xp > 0) or de- 
creasing (xp < 0) with time. 

The dynamics of the photoassociation process was 
studied by numerical solution of the time-dependent 
Schrodinger equation describing the motion of wave pack- 
ets in the ground state potential a'^S+(6s -I- 6s) (Iground) 
and in the 0~(6s ^5*1/2 + 6p ^^3/2) excited potential 
{Voyic) coupled by the electromagnetic field. The ini- 
tial state is represented by a stationary continuum wave 
function describing the relative motion (s-wave) of two 
colliding Cs (6s, F=4) atoms, with an initial kinetic en- 
ergy of 54 /iK. The wave packet in the excited potential 
T4xc is analyzed by projection onto the vibrational wave 
functions of the v' levels in the external well. The main 
results of the calculations are: 

• During the pulse, many vibrational levels are pop- 
ulated. However, after the pulse, the population 
remains only in the 15 vibrational levels resonantly 
excited around w^=98. This defines a photoassoci- 
ation window in the energy domain, for the levels 
located between ■i;^=106 and V2 = 92, bound by 
1.74 cm^^ and 3.57 cm~^, respectively. 

• In the impulsive approximation [2^ , assuming that 
the relative motion of the two atoms is frozen dur- 
ing the pulse, it is possible to analyze the numer- 
ical results by a two-state model at each internu- 
clear distance R. While the frequency of the laser 
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2.656 cm"^ 


120 kW cm"^ 


0.98 cm"^ 


15 ps 


34.8 ps 


-0 025 rm-^ns"^ 


21.2 ps 


2.32 


3.2 


v+ 


2.656 cm"^ 


120 kW cm- 2 


0.98 cm"^ 


15 ps 


34.8 ps 


+0.025 cm-ips-i 


21.2ps 


2.32 


3.5 


Vo 


2.656 cm-^ 


120 kW cm- 2 


0.98 cm"^ 


15 ps 


15 ps 





21.2 ps 


1 


0.4 


-pl22 


0.675cm-^ 


120 kW cm- 2 


0.26 cm"^ 


57.5 ps 


110 ps 


-0.002 cm-^ ps-^ 


81.3 ps 


1.91 


20 



TABLE I: Parameters for the main pump pulse V- considered in Refs. |20l l2ll : detuning maximum intensity II, energy 
range associated to the spectral width hScu, temporal widths (FWHMI) tl and rc, linear chirp parameter xpj new parameters 
used in the present work are the temporal width (FWHME) before chirping r^, and stretching factor fp (see Eq. Q in text). 
Also indicated are other pump pulses with different chirp or different detuning. In the last column, we indicate the population 
transferred to the excited state by the photoassociation pulse. 



is swept from ujl^Xp'''c to iLJL + XP'''c, the crossing 
point Rc{t) of the two field-dressed potential curves 
is moving in the [i?mirn -Rmax] range. Assuming that 
the resonance condition corresponds to a vertical 
transition at i? = Rc{t), this defines a photoasso- 
ciation window in the R domain: in the example 
chosen, i?min — 85 and i?max =107 ao. Within this 
range of distances, all the population of the initial 
continuum level can be adiabatically transferred to 
the excited state. 

In the excited state, the levels v' which are popu- 
lated have a classical outer turning point i?^^^ in 
the [i?min: ^max] range, and the two-state model 
predicts that the population transfer occurs around 
the time ty' when R^ut— Rc{tv')- For a negative 
chirp, the position of the crossing point varies from 
R max to i^mim SO that the Upper levels are popu- 
lated first (see Fig. For the pump pulse consid- 
ered, the level v'i= 106 is populated at the begin- 
ning of the time window (ti ~ (tp — rp), the level 
V2 — 92 is populated at the end {t2^ tp + tc), i.e. 
70 ps later. It is the opposite for a positive chirp, 
choosing the pulse for photoassociation. 

The total probability of population transfer to 
the excited state is reported in the last column 
of Table R] The absolute value is discussed in 
Refs. |20ll2lj|. the important result for the present 
paper being that the population obtained with Vo 
is increased by a factor 8 when a chirped pulse V- 
is considered, and by a factor of 6 when the pulse 
V}^^, with a smaller detuning, is used. 

After the end of the pulse, the vibrational wave 
packet 'i>cxc{R,t) in the excited state is mov- 
ing from the outer classical turning point toward 
shorter distances (see Fig^. We analyze it through 
partial wave packets, with main components cor- 
responding to the various levels v' . Each one is 
reaching the inner turning point at a time which 
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FIG. 1: (Color online) Photoassociation with a chirped pulse: 
For negative chirp, the level v'l is populated at the beginning 
(tl), while v'2 is populated at the end (t2) of the time window 
characterizing the photoassociation pulse. The magnitude of 
the chirp can be chosen such that, after half a vibration period 
Tyit, the partial wave packets arrive at the potential barrier 
at the same time t^. This time (ts) is optimal for radiative 
stabilization, with a second pulse, toward bound ground state 
levels. 



depends upon the excitation time t^i and the vi- 
brational period Tmb{v'). In the photoassociation 
window considered here, the vibrational periods are 
in the range 350-196 ps, yielding |Tvib(fi) — 175 ps 
and ^TyihW2) — 98 ps, respectively. For a negative 
chirp, we may choose the parameters of the pulse 
such as to create partial vibrational wave pack- 
ets which reach the inner turning point approxi- 
mately at the same time ^3 w ti + \/2T^i\,(v'^) k, 
t2 + l/2Tvib(w2)- Foi' ^ong range molecules, the vi- 
brational period is known to obey a scaling law |26j | , 
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and it is easy to choose the chirp parameter such 
that all the partial wave packets in the intermediate 
levels v'2 < v' < v[ satisfy j2ll | 



tv'iXf 



TvibK) 



tp 



(3) 



Focussing of the wave packet is then obtained at 
time ^3, after half a vibrational period. In contrast, 
for a positive chirp the lower levels are populated 
first, t2 < ti, and a spreading of the vibrational 
wave packet is expected. 

• Several authors in the past have suggested ways of 
designing laser fields for the generation of spatially 
squeezed molecular wave packets [S^, H^. Since 
they considered low vibrational levels, a numerical 
optimization procedure appeared to be necessary. 
Here, we are dealing with long range molecules, 
where the motion is controlled by the asymptotic 
potential, and the availability of scaling laws makes 
the optimization easier, as for atomic Rydberg wave 
packets 29]. 

• We should note that in the chosen example, the 
vibrational progression and the scaling law of the 
vibrational period are slightly perturbed due to 
a tunneling effect for the levels in the vicinity of 
v'=96 (23, This leads to a focussing time 
t3=tp+135 ps, slightly different from ^T^ih{v' — 
98) =125 ps. Moreover, the early population of lev- 
els which are perturbed by tunneling leads to popu- 
lation transfer to the inner well of the 0~(6s-|-6p3/2) 
excited potential. This tunneling effect is not a 
general phenomenon, but very specific to the ce- 
sium 0~ molecular curve, and to the value chosen 
for the detuning. A detailed investigation would re- 
quire consideration of two coupled channels in the 
region of the inner potential [33 , beyond the scope 
of this work. 

In the present calculations, we start from the wave 
packet in the excited state and investigate how the pop- 
ulation can efficiently be transferred to vibrational levels 
of the ground state by the second pulse. As in our pre- 
vious work, this is achieved by studying the wave packet 
dynamics by numerical solution of the time-dependent 
Schrodinger equation (cf. Ref. [23, |^ for an overview 
over numerical methods and Appendix ^ for a discus- 
sion of the representation of continuum states). 



III. THE TWO-COLOR SCHEME 

The scheme we suggest is illustrated in Fig. |21 (left), 
where we have represented, at different times, the wave 
packet created by the negatively chirped picosecond pulse 
V- ■ We consider a vibrational wave packet 'i'cxc {R, t = 
tp +tdYn), which has evolved on the excited state poten- 
tial for time tdyn after the maximum of the pump pulse 



tp. In order to discuss directly transfer probabilities, 
we will renormalize the wave packet before applying the 
dump pulse, such that at any time < >=1. 
The variation of the photoassociation probability as re- 
ported in the last column of Tabled will be reconsidered 
in the conclusion. 



A. The Franck-Condon overlap and its 
time-dependence: choice of the target vibrational 
level in the stabilization step 



We first compute the Franck-Condon overlap matrix 
elements \{'^ cxc{tdyn)\ fy>/)\'^ of the excited state wave 
packet with the bound ground state eigenfunctions (cf. 
Fig. m right). The results markedly depend upon idyn, 
and so does the probability to form stable molecules. For 
instance, just after the photoassociation pulse, at time 
tp+ 40 ps ^ tp + Tc (circles in Fig. [21 right), the Franck- 
Condon overlap factors are negligible for most bound lev- 
els of the ground state, except for the last level v" = 53 
which is very loosely bound (5 x 10^^ cm~^ US)- In con- 
trast, half a vibrational period later, at time ip-|-135 ps 
when the wave packet has moved to the inner turning 
point (diamonds in Fig. |21 right), a significant overlap is 
obtained with various deeply bound levels of the ground 
state, in particular with v" = 6, 14, 20. The large overlap 
with </'?"=i4(-^) can be explained by the coincidence of 
the inner classical turning point of the level v'—92 and 
the outer turning point of w"=14, such that w'=92 — + 
v"=14: is a Franck-Condon transition. This large value 
is restricted to a short (~ ±2 ps) time window around 
tp + 135 ps as illustrated in Fig. El where the time- 
dependence of the Franck-Condon overlap is shown. Note 
that the overlap with v" = 6 is due to population in the 
inner well of V^xc due to tunneling [33|. From Fig.O we 
can indeed see that this population remains trapped at 
the considered timescale. 

The above mentioned large overlap with the last vibra- 
tional level of the ground state just after the photoasso- 
ciation pulse {t — tp < Tc) is seen to reappear after one 
vibrational period, when the wave packet reaches again 
the outer turning point. However, we will not choose 
one of the last bound levels of the ground state as our 
target vibrational level v", since a significant population 
is already achieved in a one-color experiment [20l l2l| . 
making the implementation of a two-color scheme super- 
fluous. In fact, the specific advantage of a pump-dump 
experiment is the possibility of stabilization into deeply 
bound levels in the ground state. We will therefore fo- 
cus on populating the level v" = 14 with binding energy 
-Eb = 113 cm~^, and we will consider a second (or dump) 
pulse 2? with parameters adapted to efficient population 
transfer toward v" = 14, in particular with a time delay 
close to 135 ps. 
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internuclear distance [ ] energy [ cm ' ] 

FIG. 2: (Color online) left: The two-color scheme: A chirped picosecond pulse V- (see Table HJ excites a wave packet in the 
electronically excited state. A second, shorter pulse V (see Tables |n] llVt dumps this wave packet to the electronic ground 
state, right: Franck-Condon overlap between the time-dependent wave packet (population normalized to one, see text) in the 
excited state, and ground state vibrational levels at times t — tp+40 ps, just after the photoassociation pulse (blue circles) 
and a.t t = tp + 135 ps ~ tin (green diamonds), when the wave packet is focussed at the inner barrier; is the time when the 
photoassociation pulse is maximum. When the wave packet is focussed at the inner turning point, there is a maximum overlap 
with the level v" = 14 of the ground state. 



B. Choice of parameters for the second pulse 

For the dump pulse V, we define a central frequency 
u}^ , a maximum intensity II, a duration (FWHME) of 
the transform-limited pulse T]j, the time of maximum 
amplitude to, and a stretching factor fu characterizing 
the strength of the chirp. The choice of the second color 
is dictated by the choice of the target levels. For v" = 14 
(cf. Fig. 13 right), we have chosen a central frequency 
= 11843.5 cm~^, which is resonant with v" — 14 
for a binding energy of the excited state wave packet 
equal to 2.1 cm^^. This binding energy, computed as 
(^'cxc(i?,t)|r + V;xc|«'cxc(i?,t)), where f is the kinetic 
energy operator, stays time- independent once the pho- 
toassociation pulse is over. It is found 0.5 cm~^ smaller 
than the binding energy of the level v' = v'j^ = 98 which 
is resonant with the central frequency of the pump pulse. 
This difference reflects the non-symmetric population of 
the v' levels relative to w^, the photoassociation process 
being much more efficient for the upper levels. 

The spectral width of the second pulse V is limited by 
the requirement that mostly bound levels should be pop- 
ulated in the ground state. Since 'i;"=14 is bound by 113 
cm^^, the energy spread of the vibrational wave packet 
in the excited state can be neglected and the maximum 
spectral width is estimated as 2x113 cm""'^. Therefore, 
for short pulses V (FWHME < 93 fs), population 
transfer to the continuum is to be expected leading to 
the formation of pairs of "hot" atoms. This conclusion 



relies on overlap integrals and needs to be corrected by 
the possible energy dependence of the transition dipole 
moment /i(i?). Finally, if we wish to selectively excite 
the level v" = 14, the spectral width is limited by the vi- 
brational energy splitting around v" = 14 which is about 
7 cm~^. The typical pulse temporal width of the 
dump pulse should therefore be larger than 2 ps. In the 
following, we will discuss two schemes to form ground 
state molecules: 

1. For short pulses in the 10-300 fs range, impor- 
tant population transfer to several bound levels of 
the ground electronic state can be obtained (cf. 
Sec liVII . The goal is then to maximize the total 
population transfer. The intensity, duration and 
possibly stretching factor, i.e. chirp, of the V pulses 
will be systematically varied. 

2. For longer pulses in the picosecond domain, transfer 
of population to a single vibrational level, namely 
■i;"=14, is possible (see Sec|3). The goal is then to 
optimize selective population transfer to a chosen 
level. 

The third parameter is the time-delay between the two 
pulses, defined as tu —tp. The discussion can be guided 
by the following arguments, illustrated in Fig. ^ 

• First, we can use a simple model for the choice 
of the time delay to — tp. Considering the wave 
packet ^'cxc(^, ^d), it is possible to estimate an 
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40,3 

FIG. 3: A two-dimensional plot of the Franck-Condon overlap 
of Fig. ISIb) between the time-dependent wave packet in the 
excited state and all bound vibrational levels of the ground 
state, tp corresponds to i = 0. It is clear that the over- 
lap with deeply bound ground state levels (around v" = 14) 
is maximum when focussing of the excited state wave packet 
occurs, 135 ps or half a vibrational period after tp. In con- 
trast, when the wave packet is at its outer turning point just 
after the PA pulse and again at about tp-|-250 ps (one vi- 
brational period T^ib), the overlap with only the last bound 
levels is significant. Moreover, tunneling to and trapping of 
population in the inner potential well lead to overlap with 
ground state levels around v" = 6: this is specific to the ce- 
sium Og {P3/2) state and to the value chosen for the detuning. 
Note the stability of this population during half a vibrational 
period. 



upper limit to the probability of forming molecules 
from the sum over the Franck-Condon overlap ma- 
trix elements |(^'cxc(^, iD)|'/5^„) p. After a short 
and weak dump pulse at time to-, the wave packet 
created in the ground state is estimated by 



\E{t-tD) 



(4) 



where denotes jointly the sum over bound lev- 
els with eigenenergies E — Ei and the integral over 
continuum states with energy E, while e{E = hu) 
is the Fourier transform of the field e(t) and fx 
the dipole moment, assumed to be i?— independent. 
Eq. I0J implies that vibrational levels and contin- 
uum states which are resonant within the band- 
width e(a^) of the pulse are excited. The popula- 
tion transfer is proportional to the time-dependent 
Franck-Condon factors discussed above. Assum- 
ing a TT-pulse (i.e. a pulse which leads to popula- 
tion inversion) with bandwidth covering all ground 
state vibrational levels, the maximum probability 
is given by |(V'g(tD)|V's(^r>))P: restricting the sum 
in Eq. Q to the bound levels. The highest prob- 
abihty, 0.2, is expected for a delay of 135 ps (cf. 
Fig. when the wave packet ^'cxc(-R, + 135) is 
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FIG. 4: (top) The photoassociation pulses 'P±,Vo. (bottom) 
The upper limit to the probability of forming bound molecules 
is given by the sum over Franck-Condon overlap matrix ele- 
ments between the time-dependent excited state wave packet 
(assumed normalized to one) and all the bound levels of the 
ground state. The highest probability of about 20% is ex- 
pected for the focussed wave packet (at time t = tp + 135 ps). 



focussed after half a vibrational period. Note the 
relatively short time window in which the probabil- 
ity is maximum. A secondary maximum is visible 
after one vibrational period (250 to 300 ps) when 
the last levels are populated: much less popula- 
tion can then be transferred. Note also the revival 
at ti) ~ 400 ps, after one and a half vibrational 
period, with a larger width. The advantage of fo- 
cussing is further manifested by comparing to the 
Franck-Condon overlaps of wave packets created by 
the pump pulses V+ and Vq. The maximum prob- 
ability is now less than 10%, instead of 20% for 
the focussed wave packet, and it is spread during a 
larger time window. 



• Next, we report time-dependent calculations for 
various dump pulses where the field is treated non- 
perturbatively and we start from a focussed wave 
packet. Intensity and spectral width are varied (see 
below Sec. llVBjl in order to check the validity limits 
of the above mentioned simple model which relies 
on the impulsive approximation and weak fields. 
Then a chirp of the dump pulse is introduced to ob- 
tain a robust scheme for population inversion (cf. 
Sec. IIV C|l . The dependence on the time delay is 
studied in Sec. IIV L)l We furthermore present calcu- 
lations starting from vibrational wave packets cre- 
ated with other photoassociation pump pulses such 
(Sec. 



as V+, Vo (S ec. HVEll 
detuning (Sec. UVTT) . 



and considering a smaller 
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IV. STABILIZATION WITH A BROADBAND 
PULSE POPULATING MANY VIBRATIONAL 
LEVELS OF THE GROUND STATE 

A. Concept of a molecular tt pulse 

When the aim is to transfer as much population to 
the ground state as possible, the discussion of the pulse 
parameters is inspired by the arguments on population 
inversion in molecular electronic states developed in Refs. 
[22. .31] . The basic idea is that for pulses which are short 
on the timescale of the vibrational dynamics, the wave 
packet does not move during the pulse (impulsive approx- 
imation). The molecular system can then be reduced to 
an effective two-level system, and it is possible to find 
a TT-pulse which leads to total population transfer. A tt- 
pulse is defined as a pulse (assumed here to be transform- 
limited and with Gaussian envelope) whose Rabi angle, 

f+°° , , , f+°° -V2h[2t^ , , , 
v — iie(t)dt — / /xeo exp -j^ at (5) 

J —00 J —00 ^ 

D 

is equal to an odd multiple of tt. A pulse with Rabi 
angle equal to an even multiple of tt induces zero pop- 
ulation transfer. Introducing a chirp enforces adiabatic 
following conditions for all Rabi angles > tt and there- 
fore suppresses the oscillations in population transfer as 
a function of the Rabi angle [2ll2ll3l| . 



B. Inducing mELximum population inversion 

The time-dependent Schrodinger equation describing 
the evolution of the wave packets 'I'cxc(-R, t) and ^'g(i?, t) 
in the excited and in the ground state, coupled by the 
radiative field of the dump pulse is solved numerically. 
Our treatment of the field is non-perturbative (cf. Refs. 
[20I l2l| ') and for intensities corresponding to Rabi an- 
gles d > TT deviations from the estimate of Eq. @ are 
expected. Deviations will also occur if the impulsive ap- 
proximation ceases to be valid, i.e. for longer pulses. In 
this case, nuclear motion during the pulse cannot be ne- 
glected and the wave packet is not simply projected to 
the electronic ground state. 

For all calculations presented in this section, the pa- 
rameters of the short dump pulses are recalled in Ta- 
ble The initial condition is given by the focussed ex- 
cited state wave packet at t = tp + 135 ps, normalized to 
one such that population corresponds directly to proba- 
bility. The propagation is carried out in the time interval 
[to = tD- 4cr, tfinai =tD + 4ct] with CT = foT^ /2V2\n2. 
The final population in the ground state Pg(tfinai) is given 
by Pg(tfinai) = |(*g(ifinai)|*g(ifinai))P whilc the proba- 
bility to form bound molecules corresponds to = 

EJ(^?l*g(ifinal))P. 




Rabi angle [ units of 71 ] 

FIG. 5: (Color online) The final ground state population (top) 
and the final population of bound levels in the ground state 
(bottom) after a transform- limited dump pulse (see Table HTIl 
as a function of the Rabi angle, i.e. intensity. For all pulse du- 
rations shown, we find pronounced Rabi oscillations, i.e. the 
approximation of an efi'ective two-level systems holds. Very 
short pulses (10 fs) are able to induce complete population 
inversion, but the population which ends up in bound levels 
is limited to about 20%; the dump pulse is then dissociating 
80% of the photoassociated molecules. Longer pulses (300 fs) 
populate only bound levels and no continuum states, since, 
within their bandwidth, only bound levels are on resonance. 



In Fig. [SI we analyze the efficiency of the population 
transfer by varying the intensity of a transform-limited 
dump pulse for several durations t^, in the range of 10- 
300 fs. The approximation of an effective two-level sys- 
tem is shown to hold for all pulse durations, since oscil- 
lations as a function of the Rabi angle are clearly man- 
ifested in the final state population Pg(tfinai) and in the 
population of bound levels -Pg°""'^(tfinai)- The contrast of 
these oscillations is maximum for r^ = 10 fs, where total 
population inversion is obtained with Pg(tfinai) oscillat- 
ing between and 1. For longer pulses {t^= 300 fs), 
only 10% of the population is transferred to the ground 
state. However, an analysis of the wave packet in the 
ground state shows that for the shorter pulses, the pop- 
ulation is spread between continuum and bound levels. 
The probability to form bound molecules is at most 20% 
(t^=10 fs), so that 80% of the photoassociated molecules 
are dissociated by the dump pulse, giving pairs of hot 
atoms which usually leave the trap. With longer pulses, 
due to their narrower frequency bandwith, the transfer 
to continuum levels is no longer possible: for r^=300 
fs the 11% population transferred to the ground state 
is totally stabilized into bound levels. The maximum 
population transfer and its efficiency are reported in Ta- 
ble demonstrating that in the intermediate range of 
pulse widths, up to 100 fs, a stabilization of ~ 20% of 
the photoassociated molecules is to be expected. 

The possibility of stabilizing 20% of the photoasso- 
ciated molecules is very appealing. However, the vi- 
brational distribution of the population transferred to 
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TABLE II: Parameters for the short dump pulses considered in Sec. II VI central frequency uil^, temporal width (FWHME), 
stretching factor fo , linear chirp parameter xd , delay between the maxima of dump and pump pulse to—tp. The corresponding 
figures are reported in the last column. 



T° [is] 10 50 100 150 300 

0.984 0.746 0.374 0.195 0.114 

pbound Q 0.208 0.180 0.158 0.114 
pbound^p^ 0.216 0.273 0.487 0.811 0.999 

TABLE III: Ground state populations after a dump yr-pulse 
starting from a unity normalized wave packet in the excited 
state; ratio of population in bound levels to overall ground 
state population after 7r-pulse (cf. Fig. |^ 
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FIG. 6: (Color online) The vibrational distribution of ground 
state wave packets after vr-pulses (cf. Fig. |KJ : The distribu- 
tions after a 10 fs and 100 fs pulse are almost identical to 
the Franck-Condon overlap of the focussed excited state wave 
packet with ground state bound levels. The narrower band- 
width of a 300 fs pulse allows for less levels to be populated 
(cut-off around v" — 10 and v"=25); furthermore, no contin- 
uum states are excited (cf. Table ITTTll . 



the ground state should also be discussed, since stable 
molecules should be formed in low vibrational levels. The 
wave packet created in the ground state after a tt pulse 
is analyzed through projection on the various ip^„ (R) vi- 
brational wave functions, cf. Fig. El For short enough 
pulses (r^ =10 fs, 100 fs) the distribution follows the 
variation of the Franck-Condon overlap discussed above 
(cf. Fig. 13 right). The difference between the 10 fs and 
100 fs pulses affects mostly the population of continuum 
states (not shown in the figure). Due to the narrower 
bandwith, the longer pulse with =300 fs is only pop- 
ulating levels up to w"=25, and no continuum levels. In 
all cases the levels around v"=lA are preferentially pop- 
ulated. The efficiency of the stabilization process would 



of course be reduced by considering Rabi angles differ- 
ent from TT. Therefore, for a given pulse duration, the 
scheme is very sensitive to the intensity of the dump 
pulse. This can be avoided by introducing chirped dump 
pulses milllS- 



C. Chirping the dump pulse to suppress 
oscillations as a function of the Rabi angle 

In Fig. 13 we report numerical results obtained after 
chirping the dump pulse, starting from transform-limited 
pulses of 100 fs and 300 fs, respectively. The streching 
factor Jd is varied from 1 to 10, and the linear chirp pa- 
rameter Xd is chosen to be negative. Note that the Rabi 
angles in Fig. [T] only characterize the pulse intensities, 
they are those of the corresponding unchirped pulses (the 
Rabi angle becomes complex for chirped pulses). Fig. [7| 
clearly shows that chirping suppresses the variation of 
the final population as a function of the Rabi angle. For 
= 100 fs and /d > 5 complete population inversion 
is reached at high intensities due to non-perturbative ef- 
fects such as power broadening (top left). However, as for 
the unchirped pulses, the population of bound levels of 
the ground state remains Hmited to about 20% (bottom 
left). The advantage of chirping consists in this value of 
20% being obtained for a wide range of Rabi angles (tt 
to 4 TT in our calculations). For t^= 300 fs, nearly com- 
plete population inversion is achieved for stretching fac- 
tors fD = '2,5, and the probability of populating bound 
levels becomes close to 20%, as for the shorter pulses. 
The chirp parameter should, however, not be too large: 
for a strongly chirped pulse {Jd = 10) with correspond- 
ing transform- limited FWHM r^= 300 fs, the popula- 
tion transfer decreases, and oscillations reappear. This 
is due to the breakdown of the effective two-level sys- 
tem approximation: With a pulse stretched to 3 ps, it is 
no longer justified to neglect the relative motion of the 
nuclei during the pulse. 

In Fig. |S1 we compare the final ground state pop- 
ulation after chirped dump pulses with = 300 for 
positive and negative chirp. The stretching factor is 
again varied from 1 to 10. For pulses with correspond- 
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FIG. 7: (Color online) Overall ground state population (top) 
and population of bound ground state levels (bottom) after 
negatively chirped dump pulses as a function of the Rabi angle 
for different values of the stretch factor fo- The duration of 
the transform- limited pulses is r^= 100 fs (left) and = 
300 fs (right). Results for an unchirped pulse are plotted for 
comparison (black circles). 



1° = 300 fs 




Rabi angle [ units of 71 ] 

FIG. 8: (Color online) Overall ground state population 
(filled symbols) and population of bound ground state lev- 
els (open symbols) after negatively (top) and positively (bot- 
tom) chirped dump pulses with transform-limited FWIfM of 
300 fs as a function of the Rabi angle for different values of 
the stretch factor fo- The results for an unchirped pulse are 
plotted for comparison (black circles). 



ing transform-limited FWHM r = 100 fs, no difference 
between negative and positive chirp was observed (data 
not shown). For r^=300 fs, pulses with negative chirp 
perform slightly better than those with positive chirp. At 
high intensities (i? > 4tt), differences between open and 
filled symbols are seen, i.e. even fairly narrow-bandwidth 
pulses start to excite continuum levels. 

The effect of the chirp of the dump pulse on the final 
vibrational distribution in the ground state is analyzed 
in Fig. 1^ A chirped pulse with intensity corresponding 
to -d = n (red squares in Fig. |3 leads to very a similar 
vibrational distribution as an unchirped 7r-pulse (black 
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FIG. 9: (Color online) The vibrational distribution of the fi- 
nal ground state wave packet after moderately chirped pulses 
(/d = 5, cf. Fig. □ with T° = 100 fs (top) and = 300 fs 
(bottom) for two different intensities corresponding to 9 — n 
(red squares) and 9 — 3n (blue triangles). The vibrational 
distribution after an unchirped 7r-pulse is shown for compar- 
ison (black circles). 

circles) for both pulse durations. For higher intensities 
{•& — Stt, blue triangles), non-perturbative effects start to 
play a role, and the vibrational distribution is changed: 
For example, the maximum is not anymore observed for 
the level with the largest Franck-Condon factor {v" = 
14). Furthermore for = 300 fs, the approximation of 
an effective two-level system starts to break down due 
to the long pulse duration. While at higher intensity 
the dynamics becomes harder to interpret, the overall 
probability to form bound molecules is increased from 
18% (i9 = tt) to 21% (z9 = Btt) for t° = 100 fs and from 
11% to 16% for = 300 fs. 



D. Dependence on the time spent on the excited 
state 

After excitation by the photoassociation pulse at long 
range, the wave packet is accelerated toward shorter dis- 
tances until it reaches the inner barrier of the exter- 
nal potential well (solid blue line in Fig. ^1 top panel). 
While a small part of the wave packet tunnels into the 
inner potential well, most of it is reflected at the bar- 
rier and moves back toward larger distances (dashed 
red line). The Franck-Condon factors reflect this time- 
dependence of the excited state wave packet (Fig.^1 bot- 
tom panel, see also Fig.O. Short and sufficiently intense 
pulses (t^ = 100 fs, 1? = tt) allow for inducing maxi- 
mum population inversion in both cases, and the final 
ground state vibrational distributions completely mimic 
the Franck-Condon factors (data not shown). How- 
ever, only the wave packet at shorter distances has good 
Franck-Condon overlap with bound ground state levels, 
and ground state molecules are formed by the stabiliza- 
tion pulse. For the excited state wave packet at later 
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FIG. 10: (Color online) top: The excited state wave packet for 
tviro different time delays between pump and dump pulses: t — 
tp + Tvib/2 (solid blue line) and t = tp + T^ih (dashed red line) 
{tp denotes the time where the pump pulse has its maximum), 
bottom: The Franck-Condon overlap of the excited state wave 
packets with the ground state bound (left) and the continuum 
(right) levels. 



times, only the last two bound levels of the ground state 
are populated (except for tunneling), and a considerable 
amount of population is transferred to the continuum, i.e. 
two free atoms are created. Fig. 1101 confirms our initial 
assumption that stabilization to ground state molecules 
is optimal when the excited state wave packet is close to 
the inner turning point of the outer potential well. 



E. Dependence on the sign of the chirp of the 
photoassociation pulse 



Excited state wave packets after pump pulses with pos- 
itive {V+, red dashed lines) and negative chirp {V-, blue 
solid lines) are compared in Fig. ^] at the time of maxi- 
mum focussing at the inner barrier. The sign of the chirp 
determines how well the excited state wave packet can be 
focussed (cf. top panel of Fig. [TT|l . For positive chirp, the 
wave packet is less focussed, and the maximum probabil- 
ity to form bound molecules is 6% compared to 21% for 
negative chirp. Furthermore, almost half of these 6% are 
due to the specific topology of the 0^ -potential of Cs2, 
i.e. tunneling to the inner well. Focussing as obtained 
with V- on the other hand represents a general scheme. 
As in the previous sections the vibrational distributions 
after the dump pulse (data not shown) completely mimic 
the Franck-Condon overlap factors for pulses which in- 
duce maximum population transfer (e.g. = 100 fs, 
I? = 7r). The population of bound levels is then 4.9% for 
V+ to 18% for V-. 



FIG. 11: (Color online) top: The excited state wave packet 
a,t t = tp + 135 ps for positive and negative sign of the pump 
pulse chirp {V- and P+). (bottom) The Franck-Condon over- 
lap of the excited state wave packet with the ground state 
bound (left) and the continuum (right) levels. 
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FIG. 12: (Color online) top: The excited state wave packet 
at t = ip + Tvib/2 for two different detunings of the pump 
pulse, bottom: The Franck-Condon overlap of the excited 
state wave packet with the ground state bound (left) and the 
continuum (right) levels. 



F. Dependence on the detuning of the 
photoassociation pulse from the atomic line 

Fig- El compares the excited state wave packet for two 
different detunings of the pump pulse (V- and V^^, 
respectively). For smaller detuning (T'l^^), higher vi- 
brational levels of the excited state potential get pop- 
ulated. Due to the larger spread of vibrational peri- 
ods, it is harder to fulfil the condition, Eq. jSJ, and 
the wave packet is less focussed (cf. top panel Fig. [T^ . 
Furthermore, no tunneling to the inner well of the 0~ 
state occurs for smaller detuning. Hence, the peak in the 
Franck-Condon overlap factors around v" = 6 is missing. 
Pulses which achieve maximum population inversion to 
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FIG. 13: (Color online) left: Overall ground state population 
after transform-limited pulses with = 1 ps (black circles, 
2 ps (blue diamonds), and 3 ps (red triangles) as a function 
of the Rabi angle, right: Vibrational distributions of the final 
ground state wave packet corresponding to the points indi- 
cated by arrows in the left panel. 

the ground state can still be found (e.g. = 100 fs, 
I? = tt). However, the probability of populating bound 
ground state levels is at 3.4% considerably smaller. The 
efficiency of the complete two-color scheme is of course 
given by the efficiencies of pump and dump steps. One 
has to keep in mind that since the probability of the pho- 
toassociation step is higher for T-"!^^ (see Tabled, the 
two effects compensate and the final number of bound 
molecules is the same. 



V. HOW TO CREATE GROUND STATE 
MOLECULES IN A SINGLE LEVEL 

In order to achieve selective population of the level 
v" ~ 14, we consider longer pulses, with parameters re- 
ported in Table IIVI As explained in Sec. IIIIBI pulses 
with corresponding transform-limited FWHM > 2 
ps are expected to populate a single vibrational level. 
Fig. El reports the final ground state population af- 
ter dump pulses with = 1 ps (black circles), 2 ps 
(blue diamonds) and 3 ps (red triangles) as a function 
of the Rabi angle, i.e. pulse intensity on the left, and 
corresponding vibrational distributions on the right. A 
pulse with = 3 ps and moderate intensity indeed ex- 
cites exclusively a single vibrational level (filled red tri- 
angles in Fig. El right). The population of this level 
is Py"=i4 — 0.12 compared to the total population of 
pbound _ Q A broadening of the vibrational distri- 
bution is only observed at fairly high intensity (t? 67r, 
empty red triangles). For = 1 ps, broadening of the 
vibrational distribution sets in already at lower inten- 
sities (■)? ~ 37r, empty black circles). For — 1 ps, 
oscillations of the final population with the Rabi angle 
still persist, albeit with much lower contrast than for the 
broad-bandwidth pulses of Sec. IIVI For longer pulses, the 
impulsive approximation completely ceases to be valid, 
and the behavior cannot anymore be understood (and 
employed for control) in terms of an effective two-level 




1 2 3 4 12 14 16 
Rabi angle [ units of Ji ] level v" 



FIG. 14: (Color online) left: Overall ground state population 
after negatively (filled symbols) and positively (open symbols) 
chirped pulses with transform-limited FWHM of 3 ps as a 
function of the Rabi angle. The stretch factor is 1.5 (red 
diamonds) and 2 (blue triangles) , respectively. The results for 
an unchirped pulse are shown for comparison (black circles), 
right: Vibrational distributions of the final ground state wave 
packet corresponding to the maxima of the final ground state 
population on the left-hand side. 

system. Chirping prolongs the pulses even further, and 
its effect on the final populations and vibrational distri- 
butions is investigated in Fig. 1141 which compares mod- 
erately chirped pulses (stretching factors fo = 1.5, red 
diamonds and /d = 2, blue triangles) to the unchirped 
case (black circles) for = 3 ps. The final ground 
state population (blue triangles, fo = 2) is decreased af- 
ter a chirped pulse as compared to the unchirped case. 
However, a slight negative chirp (filled red dianidonds, 
fo = 1.5) does not significantly infiuence the population 
transfer. In general, a range of intensities corresponding 
to Rabi angles 27t < < 3tt exists where 10% or more 
of the excited state wave packet are transferred to the 
ground state, almost exclusively to the vibrational level 
v" — 14 (cf. Fig. El right). This result implies that some 
chirping (which might always occur in an experimental 
setup) does not perturb stabilization into a single level, 
i.e. the scheme is expected to be fairly robust. 

VI. CONCLUSIONS 

The formation of ultracold molecules is presently a 
very active field. Previous work has demonstrated the 
advantage of using chirped laser pulses to create ul- 
tracold molecules by photoassociation in an assembly 
of cold atoms [20l |2]| . The molecules are created in 
an electronically excited state and are short-lived. A 
stabilization step must therefore be implemented. The 
present theoretical work has discussed a two-color pump- 
dump experiment where a second pulse, delayed from the 
photoassociation (pump) pulse, transfers population to 
deeply bound levels of the ground state, creating stable 
molecules. 

Following Ref. "2^, the pump pulse (with picosecond 
duration and maximum at time tp) photoassociates Cs 
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TABLE IV: Parameters for the "long" dump pulses considered in Sec. |V| central frequency uj'j^ , temporal width (FWHME) 
streching factor fo, linear chirp parameter xo, delay between the maxima of dump and pump pulses. The corresponding 
figures are reported in the last column. 



atoms populating, via an adiabatic transfer mechanism, 
about 15 vibrational levels in the external well of the 
0~(6s ^5*1/2 + 6p potential curve. A vibrational 

wave packet is created which is bound by 2.1 cm^^ and 
localized in the region of the outer turning point. After 
the pulse, the wave packet moves toward the inner bar- 
rier of the potential curve (located at i?in). The chirp pa- 
rameter of the pump pulse can be chosen such that after 
half a vibrational period, at time tp + Tvib/2, this wave 
packet is focussed at the barrier '2(J|. Our analysis of 
the Franck-Condon overlap integrals has shown that the 
overlap with deeply bound levels v" of the ground triplet 
state is then maximum. A vertical transition toward the 
level v"=14:, bound by 113 cni^^, is even possible, since 
the outer turning point of the corresponding wave func- 
tion coincides with i?j„. We therefore have discussed how 
population can efficiently be transferred to the ground 
triplet state by a short dump pulse which should be in 
the femto- or picosecond domain and optimally delayed 
from the pump pulse by T^ih/'i- 

The validity of this intuitive model has been checked 
with numerical calculations. We have considered two 
kinds of pulses: broadband pulses populating many vi- 
brational levels of the ground state, with the aim of max- 
imizing the population transfer, and longer pulses for se- 
lective population of a single level, v"—14:. 

The first set of calculations has been performed for a 
time delay of Tvib/2. We have shown that for short pulses 
and a Rabi angle equal to tt, total population inver- 
sion can be obtained, i.e. the vibrational wave packet in 
the excited state is completely transferred to the ground 
state. Analysis of the wave packet shows that a max- 
imum of 20% of this population goes to bound levels. 
The remaining 80% correspond to continuum states, i.e. 
the dump pulse also creates pairs of hot atoms. The vi- 
brational distribution reflects the variation of the Franck 
Condon overlap. Oscillations in the population trans- 
fer as a function of the Rabi angle (i.e. laser intensity 
for a fixed pulse duration) can be suppressed by chirp- 
ing the dump pulse. Thus, total population inversion is 
obtained for a wide range of intensities (corresponding 
to Rabi angles of tt to Att in our calculations) with an 
effective 20% population transfer to bound levels. Em- 
ploying a chirped dump pulse therefore yields a robust 
scheme in which the probability to create stable ground 



state molecules is 20%, independent of the exact value of 
laser intensity. 

By varying the time-delay between the two pulses, a 
strong variation in the population transfer which follows 
the variation of the Franck-Condon factors has been ob- 
served. Only at shorter distance, the wave packet in 
the excited state has good Franck-Condon overlap with 
bound ground state levels, and ground state molecules 
are formed by the dump pulse. The highest efficiency is 
achieved in a narrow time window (Tvib/2 ± 2ps) during 
which the wave packet is maximally focussed at the bar- 
rier. For different time-delays, only the last two bound 
levels of the ground state are populated, and a consider- 
able amount of population is transferred to the contin- 
uum, i.e. two free atoms are created. Note that without 
a second pulse, stabilization would occur via spontaneous 
emission on a nanosecond timescale. The probability for 
creating ground state molecules is then given by a time 
average. Since the wave packet spends more time at large 
distances than close to the barrier, times with unfavor- 
able Franck-Condon overlaps dominate this average. It 
is therefore essential to employ a second pulse to fully 
take advantage of the time- dependent probability to cre- 
ate bound ground state molecules after photoassociation 
with picosecond pulses. 

The calculations have confirmed our intuition that sta- 
bilization to ground state molecules is optimal when the 
excited state wave packet is close to the inner turning 
point of the outer potential well. The best localization of 
the wave packet at the inner barrier is brought about by 
negatively chirping the pump pulse. There is yet another 
advantage when using a negative chirp of the pump pulse: 
The photoassociation probability is increased by a factor 
of 8 as compared to a transform-limited pulse. Further- 
more, the probability of creating bound molecules by the 
dump pulse is enhanced by roughly a factor of 3 for a 
negatively chirped pump pulse as compared to positively 
chirped or transform-limited ones. The number of sta- 
ble molecules formed in a pump-dump sequence can be 
therefore be increased by 24 when appropriately choosing 
the pulse parameters. 

We have also considered a pump pulse with smaller 
detuning, for which the photoassociation probability is 
larger. However, due to a decrease of the Franck-Condon 
overlap, even for the maximally focussed wave packet, 
the stabilization step is much less efficient. Therefore, 
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the total number of stable molecules is roughly the same 
as in the previous case. 

Finally, we have employed longer dump pulses with a 
duration of about 3 ps to selectively populate the u"=14 
level in the ground state. Up to 12% of the photoassoci- 
ated molecules can then be stabilized into a single level, 
which is fairly efhcient. Note that these molecules are 
created coherently by the pump-dump sequence. 

The full discussion of the definition of the optimal pulse 
should include the possibility of a repetition scheme, con- 
sidering further pump and dump pulses. In treating a 
repetition scheme, we also have to take into account how 
the pump and dump pulses modify the continuum wave 
function describing two colliding atoms. This is beyond 
the scope of the present study. 

In conclusion, a two-color pump-dump scheme is rather 
promising for the creation of stable cold molecules. Due 
to the choice of the Cs2 0~ potential for the photoasso- 
ciation step , the present work has analyzed formation of 
molecules in the i;"=14 level of the ground triplet state, 
bound by 113 cm~^. Many other schemes, e.g. [Tsllssl l. 
are possible and should be considered to bring population 
into the v"=0 level of the ground state which will open 
the way to the formation of stable molecular condensates. 
In this respect, we would like to point out the generality 
of the two-color pump-dump scheme which only requires 
efficient photoassociation and stabilization mechanisms 
invol ving either long range wells in the excited poten- 
tials lU or resonant couphng |0, 0, 0|, but is not 
restricted to a particular system. 
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APPENDIX A: A NOTE ON THE 
REPRESENTATION OF CONTINUUM STATES 
WITH A MAPPED FOURIER GRID 

The Fourier grid employs a finite number of grid points 
and periodic boundary conditions to represent a quan- 
tum wave function. For V{R) = 0, this corresponds to 
the particle in a box problem. The scattering continuum 
is hence represented by a finite number of discrete box 
states. The highest momentum which is correctly rep- 



resented on the grid is given in terms of the (constant) 
spacing AR between two grid points, Pmax — tt/AR. 

In order to model long-range interactions which require 
large grids but not small grid steps, a mapping procedure 
is introduced. The grid steps in the asymptotic region 
become very large, correspondingly high momenta are 
cut off. This alters the representation of continuum states 
on the grid. When diagonalizing the Hamiltonian, two 
kinds of continuum eigenfunctions are obtained: true box 
levels and wave functions with an unphysical shape. The 
latter will be called "cutoff-states" in the following. 

The box levels can be easily identified by the scaling 
of their eigenenergies. En ^ (except for the first few 
states just above the dissociation limit where the energies 
are perturbed by the threshold) . The corresponding wave 
functions are standing waves of the box. 

The eigenenergies of the cutoff-states approximately 
scale with some high power (~ 20) of the box index, the 
highest eigenenergy being determined by the largest al- 
lowed grid step which is a parameter of the mapping. 
The corresponding eigenfunctions show very fast oscil- 
lations and don't obey the boundary conditions at i?inin 
and i?max- It can nevertheless be shown that these cutoff- 
states constitute an effective representation of the true, 
physical continuum for energies higher than the energy of 
the last true box level. To this end, the Franck-Condon 
overlaps of an excited state wave function with the eigen- 
functions of all ground state levels (i.e. bound, box and 
cutoff-levels) are calculated and compared (i) for different 
mapping parameters and (ii) to Franck-Condon overlaps 
obtained by Numerov integration ,35]. 

In the first set of calculations, the largest allowed grid 
step, AiJmax, is varied. The Hamiltonian is diagonalized 
for each A7?„iax and an increasing number of grid points, 
Nii. For each Ai?niax, the number of cutoff-states turns 
out to be independent of Nji, i.e. increasing Nji leads to 
a larger number of true box states. In particular, for the 
smallest ARj^ax and largest Nn used (Ai?niax = 0.33 ag, 
Nfj — 8191), a continuum representation is obtained 
which contains almost no cutoff-states (168 vs. 7969 true 
box states). The excited state wave function is then in- 
terpolated to these grids, and the Franck-Condon over- 
laps are computed. Note, that the Franck-Condon over- 
laps for continuum states need to include the correct en- 
ergy renormalization The Franck-Condon overlaps 
turn out to be roughly independent of the grid param- 
eters. We are particularly interested in the sum over 
the Franck-Condon overlaps which determines the prob- 
ability of dissociating the excited state molecules into 
hot, free atoms with a short dump pulse (cf. Sec. IIV|I . 
This sum is roughly constant for all grids considered, 
varying between 0.21 for Ai?max = 33 ao and 0.17 for 
Ai?inax = 0.33 ao. The variation is caused by the in- 
terpolation of the excited state wave function to the new 
grids which is not exact, the approximation getting worse 
for decreasing AR^^ax- 

In a second set of calculations, Franck-Condon overlaps 
with all ground state levels were calculated for selected 
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excited state vibrational levels with both the mapped 
Fourier grid and the Numerov method. Oscillations like 
those in Fig. ^| are observed. The resolution of these 
oscillations in the mapped Fourier grid method depends 



on the number of grid points, Nn. Good agreement be- 
tween the two methods is obtained for Nji — 1023, the 
value used throughout Sees. II VI andlvl 
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